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ANALYSIS OF GROUND REFLECTION OF JET NOISE OBTAINED 
WITH VARIOUS MICROPHONE ARRAYS OVER AN ASPHALT SURFACE 

by J. H. Miles 
Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio 

ABSTRACT 

Ground reflection effects’ on the propagation of jet noise over an 
asphalt surface are discussed for data obtained u§ing a 33.02 cm (13-in.) 
diameter nozzle with microphones at several heights and distances from the 
nozzle axis. 

Analysis of ground reflection effects is accomplished using the con- 
cept of a reflected signal transfer function which represents the influence 
of both the reflecting surface and the atmosphere on the propagation of the 
reflected signal in a mathematical model. The mathematical model used as a 
basis for the computer program was successful in significantly reducing the 
ground reflection effects. The range of values of the single complex num- 
ber used to define the reflected signal transfer function was larger than 
expected when determined only by the asphalt surface. This may indicate 
that the atmosphere is affecting the propagation of the reflected signal 
more than the asphalt surface. Also discussed is the selective placement 
of the reinforcements and cancellations in the design of an experiment to 
minimize ground reflection effects. 
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INTRODUCTION 


Acoustic data takmg programs for experimental res ear cli ou axi.ui.aft 
noise reduction using ground facilities often have problems involving 
interpretation of data because of ground reflections. Acoustic reflections 
from the ground cause problems, since in some experiments, reflections dis- 
tort the true sound spectrum and this can lead to incorrect conclusions 
about spectral shapes and overestimation of the overall sound pressure 
level (OASPL) relative to the free field QASPL. This is particularly evi- 
dent with large-scale and full scale test rigs where, due to the source/ 
microphone geometry, the peak sound pressure level (SPL) values in the 
spectra can be in the midst of that part of the spectrum where the magnitude 
of the distortion caused by the interference of the direct and reflected 
pressure wave is a maximum and the resulting cancellations and reinforce- 
ments of the pressure wave are strongest. 

Many different microphone arrangements have been used at NASA Lewis 
acoustic test facilities. Experience shows that centerline horizontal 


microphone arrays of the order of 1.22 to 3.66 m (4 to 12 feet) above the 
ground are subject to such ground reflection problems. If ground micro- 
phones are used, subject to certain environmental conditions (such as the 
absence of temperature inversions) being met, the low frequency acoustic 
data are not subject to ground reflection reinforcements and cancellations; 
However the high frequency regions may be subject to ground reflection prob- 
lems. Overhead-boom microphone array have also been used with and without 
ground treatment to reduce ground reflection problems. The use of such a 
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boom is generally prohibitive for large or full-scale engine tests because 
of the 30,48-60.96 m (100-200 foot) microphone radius necessary to measure 
the engine nosie. 

Considerable effort has gone into experimental and theoretical inves- 
tigation aimed at minimizing reflection effects. Early work at NASA Lewis 
(ref. 1) has served as the theoretical basis for much succeeding work 
(refs. 2 to 4). The investigations of references 1 and 4 provide under- 
standing of the effect of a change in microphone geometry on the measured 
data for many idealized cases. However, experience in the study of broad- 
band jet noise at NASA Lewis has shown that the direct application of the 
results from the theoretical and experimental studies described in refer- 
ences 1 and 4 to practical non- ideal cases can cause difficulties. In 
many cases, the correction based op the studies of references 1 and 4 would 
lead to a more distorted spectral shape than the one measured. 

The present paper discusses acoustic spectral data from cold-flow noz- 
zles taken with microphone arrays at a variety of distances and heights 
above an asphalt surface, A method of analysis is then applied to these 
data to account for the ground reflection problems in the data. The method 
calculates the lossless spectra from the measured data. In addition, the 
reverse is also possible; i.e. the "rig" spectrum can be predicted from a 
free field spectrum. This makes possible selective placement of reinforce- 
ments and cancellations by using the computer method to identify desirable 
microphone positions. Necessary analytical expressions and procedures based 
on earlier work are included. 
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The paper will first present the equations used to model the jet 
broadband spectrum and the ground reflection. The calculation method 
will then be discussed. Finally a comparison of this approach with ex- 
perimental data will be presented. 

ANALYSIS 

Simulation of Measured Spectrum 

The present analysis considers a mathematical model for the measured 
broadband subsonic jet noise third octave sound pressure level (SPL) spec- 
trum generated by a convergent nozzle over an asphalt surface* 

A schematic of the nozzle microphone geometry relative to the 
asphalt is shown in figure 1. The parameters in this model can be adjusted 
so that the resulting calculated spectra is a least square fit to the mea- 
sured spectra. From the model the following are achieved: (1) corrections 

are obtained that enable the calculation of the lossless spectrum from the 
measured data and (2) an equation for a smooth least square curve fit to 
the lossless spectrum is generated. 

The equation representing the SPL spectrum are written as the sum of 
a lossless third octave spectrum level, SPL + , a ground reflectance, | T | ^ , 
and an atmospheric attenuation level, A. Thus, in dB 

SPL cal - SPI 4l + 10 lo SlO |T| 2 + A (1) 


referenced to a pressure p r where 
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p r = 2x10 Pascal (2x10 ^ y bar) 

(All symbols are defined in the Symbol List Appendix A.) 

The method of analytically representing the lossless SPL spectra, 

SPL^ai used herein was briefly described in reference 2 and presented in 
reference 3 in more detail. The necessary equations used herein from 
these references are presented in Appendix B. Briefly, the SPL aa ^ in any 
third octave band is determined from: (1) an overall sound pressure level 
value; (2) an effective diameter, D e , and an effective velocity, V e , used 
with frequency to determine the Strouhal number, St, (St = f D e /V e ) ; and 
(3) a pressure spectral density level curve, PSD*^, which is a function 
of Strouhal number and which is normalized to have unit area. Thus 

OASPL + /10 D 

SPL+ al - 10 log l0 10 P Af c + PSD+ al (B-l) 

e 

The ground reflectance, [t|^, is determined by two factors. The first 
factor is the complex reflected signal transfer function written in terms 
of its magnitude Q(s) and its phase 6(s) as 

G(s) = Q(s) eJ S(s) (2) 

where s is the Laplace transform variable (s = j27Tf). The second factor 
is the microphone and test site geometry (fig. 1). 

The model used to take account of the effect of microphone and test 
site geometry assumes that the noise source is distributed over a region 
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and represents this by taking into account signals from N independent 
sources. This model creates a spectrum dependent on a range of time 
delays. The time delays are calculated by assuming the N sources are 
distributed vertically about the centerline at the nozzle exit. This 
assumption was used to study acoustic interference of the noise produced 
by a single jet due to reflections from a surface (ref. 4) and in refer- 
ences 2 and 3 to study acoustic interference of flap noise due to reflec- 
tion from a surface. 

The equation for the ground reflectance for third octave spectra 
taken from references 2 and 3, is given by 


T(f c ) 


= 1 + 


Q 2 ( f c ) 


2 Q( f c>f 

+ N z “cos 
n=l 


[27rf c T n - 6 (f c )] 


sin7rAf i 
c n 

TrAf 

c n 


(3) 


where the functional dependence of Q and <5 has been changed from the 
Laplace transform variable s to f c to indicate that Q and 6 are 
assumed constant over a third octave band having center frequency, f c . 

In equation (3) the time delay, T , between the reflected and direct 
pressure signal is given by 

s 9 (n) - s 1 (n) 

T = . 

11 

The symbol z is the mean value of the ratio of the reflected signal 
travel distance to the direct signal travel distance given by 


s 2 (n) £r 2 + 

jh s (n) + h 0 ] 

J 

s l (n) {r 2 + 

OT 

/’“N 

0 

1 

S3* 

. ° , 

JT* 

CM 


z 


(5) 
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Empirical formulas for computing the atmospheric attenuation level 
A(f c ,r) can be found in references 5 and 6. The data reported in this 
paper were corrected using the atmospheric attenuation level calculated 
by the method of reference 6. 

Reflected Signal Transfer Function 

k 

The reflected signal transfer function, G(s), is used to represent 
the effect of reflections and it arises as follows. The reflected signal 
is assumed to be related to the direct signal by a convolution integral of 
an impulse response with the direc.t pressure signal. This impulse response 
represents the effects of both the reflecting surface and the atmosphere on 
the propagation of the reflected signal. The Laplace transform of this 
impulse response is the reflected signal transfer function, G(s), (eq. (2)). 
This approach is presented in more detail in references 2 and 3. 

The magnitude, Q(f c ), and phase <5(f c ) of the reflected signal transfer 
function, G(f c ), must be specified in order to calculate the reflectance, 

]t (f c )j^. Neglecting atmospheric effects the value of the reflected ray 
transfer function has been determined from physical assumptions concerning 
the nature of the acoustical signal i.e. plane wave or spherical wave and 
the nature of the boundary at the relevant surfaces i.e. scattering from a 
rough surface, absorption by an infinite surface, or a multilayer surface. 
However, realistic (i.e. non-plane wave) mathematical models cannot be 
solved exactly and always requires approximations. 

The empirical approximation found most useful in specifying G(f c ) was 
to assume that Q(f c ) and 6(f c ) are constant over the frequency range of 
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interest. Thus 



( 6 ) 


A discussion of this approximation follows. The reflected signal transfer 
function, G(f c ), has not been related to any physical model in the previous 
sections. However if the surface material can be characterized by a com- 
plex ground surface impedance an d a speed of sound C2» then neglect- 

ing atmospheric effects, G is the plane wave reflection coefficient, Cr, 


Z2 cos ©2 - cos ©2 
7 ft Z ft" = G = Q e 


^R cos ©^ + cos ©2 


where 


and 


sin ©^ sin ©2 



( 8 ) 


- P c Q ( 9 ) 

and the plane wave is incident at an angle ©^ to the normal of the sur- 
face (refs. 7 , 8). 

For a spherical wave the same reflection coefficient is applicable 
for source-microphone distances much larger than 3 wavelengths. This is 
shown for the analogous electromagnetic case (ref. 9 ). This result as 
applied to the acoustic problem was previously discussed in reference 10 
and used in reference 11 . Reference 11 also discussed the general spheri- 
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cal wave reflection problem using results of the analogous electromagnetic 
case obtained in reference 12. 

The ground absorption coefficient ag is related to the reflection 
coefficient by 

“g = 1 - |c R | 2 (10) 

Thus, the square of the magnitude of the reflected signal transfer 
function is related to the absorption coefficient by 

|G(f c )| 2 = | Q[ 2 = | C R | 2 = 1 - a g (11) 

i 1 2 

Hence, assuming that |G(s)| is independent of frequency corresponds to 
assuming that the absorption coefficient a is also independent of fre- 
quency. 

The author is not aware of data supporting the assumption that Q(f c ) 
and <$(f c ) are constant over the frequency range of interest for asphalt 
surfaces. However impedance tube measurements by Dickinson and Doak of 
the normal absorption coefficient of sharp sand over a frequency range 
from 200 to 1000 hertz for a moisture content ranging from 0 to 10 percent 
Indicated the normal absorption coefficient was not strongly dependent on 
frequency at a given moisture content (ref. 13). 

Although it is not known if asphalt behaves like sharp sand, the 
assumption will be made as a necessary part of the procedure that Q(f c ) 
and 6(f c ) are constant over the frequency range of interest. 
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Determination of Lossless Spectrum 


Tlie lossless spectrum is calculated by curve fitting equation (1) to 
measured data by adjusting the parameters that determine | T | ^ and 
SPL^ a ^. A least square fit to the lossless sound pressure level is then 
given by equation (B— 1 ) for SPL^^ . 

The lossless sound pressure level can also be calculated by subtract- 
ing the atmospheric attenuation level and the reflectance from the measured 
data. Thus 


SPL^; - SPL^ - 10 log 10 | T| 2 - A 


( 12 ) 


Another procedure for obtaining lossless spectra is to use ground 
microphones. For a microphone on the ground the time delays are close to 
zero since the reflected signal path differs very little from the direct 
signal path. Thus from equations (3), (4), and (5), over a large frequency 
range, 

1 T(f c ) | 2 = (1 + Q) 2 (13) 


and the reflectance is independent of frequency. The maximum value of Q 
is 1. For this value of Q, 

10 log | T (f )| 2 = 6 dB (14) 

10 c 

Thus, if the value of Q is 1, the lossless spectrum could be calculated 
from a ground microphone measured spectrum over a large frequency range by 
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SPL = SPL - 6 dB - A 
m m 


( 15 ) 


For a surface characterized by a complex impedance Z 2 much larger 
than the impedance of air (Z^ = P c Q ) equation (7) indicates that Q 
will be 1, and equation (15) can be used to calculate the lossless spec- 
trum from a ground microphone measured spectrum. 

CALCULATION METHOD 

The complete equation used to simulate the SPL at a microphone is 
found by substituting into equation (1) equations (B-l) , (B-3) , and (3). 
Thus the analytical representation is 


OASPLp/lO 

SPL cal “ 10 lo SlO 10 Af 


c V 


J^d r -i J^b r 1 

c o + E E [ u(i) > n P (i) ] - E E l u(k) > n P (k) J 

X— 1 1c— 1 


( 16 ) 


+ 10 log 


10 


2 N 

1 + 2 + N z £ C ° S (2lTf c T n ' &) 

z n=l 


sin (TTAf c T n ) 
irAf„ T 

C n 


+ A (f c ,r) 


A computer program is used to iteratively adjust the parameters 
0ASPL + » C 0 , Stp(i), rip (i) St p (k), rip(k) which determine the spectral shape 
and the parameters Q and 6 with determine the ground reflectance. A 
diagram of the method is shown in figure 2. The calculation method is 
described briefly in reference 2 and presented in more detail in reference 3. 
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The parameters are adjusted to minimize the total "cost", C T , consisting 
of the sum of the squares of the differences between measured and calcu- 
lated SPL values, C e , and a "penalty cost", Cp, (refs. 14 and 15) due to 
selection of parameters that violate the parameter constraints. (See eq. 

(B-2) for 0ASPL + constraint.) The total "cost" function is nonlinear. 

P 

The problem of minimizing the nonlinear "cost" function was solved by 
search techniques. These nonlinear minimization techniques consist of 
systematic procedures for varying the "cost" function parameters until a 
minimum value of the "cost" function is found. The basic nonlinear minimi- 
zation methods available are discussed in references 14 and 15. To provide 
flexibility in the choice of M a and and due to the presence of 

parameter constraints, a search technique which does not require evaluation 
of derivatives was chosen. The nonlinear minimization search technique was 
that of Powell (ref. 16). The nonlinear minimization computer program used 
was adapted from reference 17. The optimization procedure was previously 
applied successfully to a similar identification problem in references 18 
and 23. 

The input information for the method consists first of the constants 
required for the reflection model. These consist of the following: the 

number of sources N, the mean value of the ratio of the distance the 
reflected ray travels to the distance the direct ray travels, Z; and the 
time delay T n associated with each of the N sources. The number of 
time delays, N, was selected empirically to be 5 because for the range of 
time delay values used increasing the number of delays beyond five had no 
further effect on the third octave spectrum. 
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Second, the constants M & and for the PSD ca -^ model are needed 

(see ref. 3). For the data considered, was selected to be 2 and M a 
was selected to be 1. Last, initial values of the following parameters 
are necessary. 

OASPlf, c 0 , Stp(i) and n p (i); i = 1, St p (k) and n p (k); k - 1,2; Q and 6 

Procedures for selecting initial values are discussed in reference 3. 

The computer output is the final values of these parameters. Using 
the parameters OASPL^, C Q , St p (i) and Ti p (i); i = 1, St p (k) and r| p (k); 
k 53 1,2 in equation (B-l) and (B-3) the lossless spectrum SPL^^ can be 
computed. From the value of the atmospheric attenuation and using the 
reflectance parameters Q and <$ in equation (3) the lossless spectrum 
SPL^ can be computed from equation (12). 

APPARATUS 

Noise tests were conducted with a 33.02 cm (13 inch) diameter circular 
convergent nozzle having the nozzle centerline 3.89 m (12 3/4 feet) above 
grade. For this rig, dry cold air (at about 520° R) was supplied to a 
40.64 cm (16- inch) diameter gate valve from the Center’s 862 kPa (125 psia) 
air supply system by way of a 60.46 cm (24-inch) diameter underground pipe 
line. Details of the flow and muffler systems are given in reference 
19. 

The noise data were measured with sixteen 1.27 cm condenser micro- 
phones with windscreens located as shown in Table I, The microphones were 
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arranged to be approximately 145° from the nozzle inlet at a radial dis- 
tance of 7.62, 15.24, and 22.86 m (25, 50, and 75 feet) from the nozzle 
exit. The mikes near the surface were thus moved closer to the nozzle 
and at larger angles at the 7.62 and 15.24 m (25 and 50 feet) radius. 

The air flow rate was measured by an orifice flowmeter located in a 
straight section of the underground air supply line upstream of the gate 
shutoff valve. Pressure drop across the orifice flowmeter and static 
pressure upstream of the flowmeter were measured by strain-gage pressure 
transducers. Strain-gage pressure transducer were also used to measure 
total and static pressure upstream of the nozzle. All pressures were re- 
corded on strip-chart instruments. 

Weather data were also monitored and/or recorded (i.e. barometer, 
temperature, humidity, wind speed and direction). 

The noise data were analyzed on-line with an automated 1/ 3-octave 
band spectrum analyzer. The analyzer determined sound pressure level spec- 
tra (referenced to 2xl0~^ p„ (0.0002 microbar)) between 50 and 20,000 Hz at 
each microphone position. Three noise samples were taken at each microphone 
and treated to reject background disturbances and random errors and to ob- 
tain an average value. 

Noise measurements were made for a series of nominal nozzle pressure 
ratios. These were 1.1, 1.2, 1.3, 1.4, 1.5, 1.7, 1.8, 2.0, and 2.1. The 
static pressure in the free jet stream was assumed to be atmospheric, and 
the the total temperature was assumed to be the same as the temperature 
measured upstream of the nozzle. Jet exhaust velocities were calculated 
from the isentropic equations. 
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The test procedure was to obtain a steady flow condition for a given 
total pressure upstream of the nozzle- Noise data were then taken at each 
microphone location in succession. After three complete acoustic measure- 
ment cycles the airflow rate was changed and the process repeated. 

RESULTS AND DISCUSSION 
Sample Results 

The computer calculation procedure provides the three major quanti- 
ties that are added together as indicated by equation (1) to calculate the 
spectrum, SPL ca ^. The computer calculation method selected the parameters 
(OASPLp, C 0 , Stp(i) and U p (i); i = 1; St p (k) and r) p (k); k = 1,2; Q and 6) 
thereby adjusting equation (16) such that the shape of the SPL Qa ^ spectrum 
is in good agreement with the measured data, SPL m . 

Figure 3 shows typical results obtained by applying the computer method 

described herein to the data. This data and all following data are from the 

i 

33-02 cm (13 inch) diameter conical nozzle operating at exhaust velocity of 
301 m/sec (988 ft /sec). The microphones were- at a nominal directivity angle 
of 145° from the nozzle inlet. 

A measured spectrum is compared with a computer calculated spectra in 
figure 3(a). Figure 3(a) also shows the calculated lossless spectrum, 
SPL^ a ^, which as indicated previously is based on a least square curve fit 
to the measured data. However while it is a smooth curve (since it takes 
into consideration all the data and is based on a specific formula for the 
spectrum) , the resulting spectrum shape can have sections which do not have 
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the smooth convex curve shape conventionally associated with jet noise 
spectra. 

In order to obtain a spectrum which has an unbiased shape, it is 

+ 

useful to calculate a lossless spectrum, SPL m , by applying the calculated 

ground reflection correction and the atmospheric attenuation correction 

to the measured spectrum data as indicated by equation (12). However, the 

resulting curve is in general not smooth. Figure 3(b) compares calculated 

analytical lossless spectrum SPL^ a ^ and calculated measured lossless spec- 

trum SPL from the measured data shown in figure 3(a). 
m 

J 

A SPL curve smoothed by fairing a curve to either the SPL ^ or the 

SPL + curve of figure 3(b) is shown in figure 3(c) along with the SPL + 
m cal 

curve and SPL + curve. Fairing to the SPL + , curve may be done to remove 
m cal 

any bias in the least square curve fit due to (1) selection of the func- 
tional form equations and/or (2) the limited number of parameters used to 

curve fit with. Fairing to the SPL curve may be done to remove irregular- 

111 

ities introduced or not removed by the calculated reflectance and/or to 
remove other irregularities thought to be present in the data. 

In figure 3(d), the final lossless spectrum curve, SPL + , obtained by 
fairing is compared with the original spectrum to show the effect of the 
combination of the two corrections (atmospheric absorption and ground re- 
flection) . Also shown is the original SPL spectrum corrected by removing 
the atmospheric absorption. The SPL m corrected for atmospheric absorption 
for these data is shown to be about 1 dB above the lossless spectrum from 
5000 hertz to 10,000 hertz in accordance with correction for the ground 
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reflection. The atmospheric attenuation loss above 10,000 hertz causes 
the measured data, SPI^, to be lower than the free field SPL + curve. 

Hereinafter, for simplicity the results of the computer method (i.e. 

+ + 

the computer corrected SPL^ spectra and the equation calculated SPL ) 

will be shown rather than smooth SPL + curves produced by fairing to the 

SPl/t or SPL + curves. 
m cal 

Results for Various Microphone Heights 

Measured spectra . - Typical sound pressure level spectra for the 
33.02 cm (13-inch) diameter nozzle taken at various microphone heights at 
distances of 7.62, 15.24, and 22.86 m (25, 50, and 75 ft) from the nozzle 
exit are shown in figure 4. Cancellations and reinforcement effects due to 
ground reflections are seen to be quite pronounced for frequencies below 
around 1500 Hz. The data also show the wide range of pressure level values 
that can be obtained in this frequency range due to variation in microphone 
height. The data presented are for a nozzle exhaust velocity of 301 ra/sec 
(988 ft/sec). Similar trends were observed at other velocities. 

Lossless spectra . - Lossless SPL spectra were calculated based on 
correcting the measured values (SPL m , eq. (12)) for the data of figure 4 
for the same exhaust velocity and microphone heights and distances. For 
the data taken with microphones on the ground, two different procedures 
were used to correct the ground microphone data for reflection effects: 
the first was to subtract 6 dB; and the second was to subtract a reflec- 

2 _i_ 

tance, | T j , found by the computer method. The results for the SPL^ spec- 
tra are shown in figure 5. 
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A value for each measurement of the magnitude of the reflected signal 
transfer function Q and its phase 6 was found with the computer method 
described herein. The values of Q are listed for each calculated spec- 
tra in figure 5. The resulting range in values of Q could be attributed 
to non-uniformity in the asphalt absorption and/or atmospheric effects on 
sound propagation. The variation of Q is discussed in Appendix C. 

The results shown in figure 5 indicate that, in comparison with fig- 
ure 4, the scatter in the data has been greatly reduced in the region where 
ground effects were pronounced. At the various distances from the nozzle 
data at each of the microphone heights (continuous lines) are in relatively 
good agreement and indicate the relative effectiveness of the correction 
calculation. 

Examination of detailed variations. - Figure 5(a) shows that the SPL 

° m 

spectra for microphones 7.62 m (25 ft) from the nozzle and heights of 0.30, 
.91, 1.52, and 3.89 m (1, 3, 5, and 12 3/4 ft) are in good agreement over 
the entire frequency range. The 7.62 m (25 ft) radius ground microphone 
SPL m spectra corrected using a value of Q = 0.91 for the magnitude of 
the reflected ray (Jt|^ - 5 dB) is consistent with the SPL^ data at 0.30, 

0.91, 1.52, and 3.89 m (1, 3, 5, and 12 3/4 ft) up to 1600 hertz. No 
explanation is available, at present, as to why there is no agreement above 
1600 hertz. 

Figure 5(b) shows that the SPI/j^ spectra for microphones 15.24 m (50 ft) 
from the nozzle and at heights of 0.91, 1,52, and 3.89 m (3, 5, and 12 3/4 ft) 
are again in good agreement. However, above 5000 hertz the SPL+ spectra for 
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0.91 (3 ft) is higher than the 1.52 and 3.89 m (5 and 12 3/4 ft) data. 

The 15.24 (50 ft) ground microphone SPL+ data shown in figure 5(b) cor- 
rected by subtracting 6 dB are in agreement with the SPL^ data taken at 
0.91, 1.52, and 3.89 m (3, 5, and 12 3/4 ft) only up to 500 hertz. At 
higher frequencies than 500 hertz, subtracting 6 dB appears to produce 
too large a correction. The ground microphone data corrected using Q - 
0.637 ( j T| 2 - 4.3 dB) appears to be undercorrecting the ground microphone 
spectra below 500 hertz; however, above 500 hertz the use of the calculated 
correction causes the ground microphone data to be in agreement with the 
SPL^ spectra for microphones at 0.91, 1.52, and 3.89 m (3, 5, and 12 3/4 ft). 

Figure 5(c) shows that the SPL* spectra for microphones 22.86 m (75 ft) 
from the nozzle and heights of 0.91, 1.52, and 3.89 m (3, 5, and 12 3/4 ft) 
are in very good agreement over the entire frequency range. The ground 
microphone SPL^ data corrected using Q = 0,637 ( | T ] 2 - 4.3 dB) is in good 

agreement with the data up to 10,000 hertz. However, the 22.86 m (75 ft) 
ground microphone SPL^ data shown in figure 5(c) corrected by subtracting 
6 dB appears to be overcorrected. 

Individual comparisons . - Comparisons of the measured sound pressure 
level and the lossless sound pressure level calculated by applying the cor- 
rections found by the computer method at each of the four heights are shown 
in figure 6. The reduction of the magnitude of the effect of the reinforce- 
ments and cancellations due to ground reflection in each case is seen more 
clearly in the figure. 

As indicated previously, the method also provides a least square curve 
fit to the SPL^ level spectra. Figure 7 shows comparisons between SPL^ and 
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SPL^ al for two cases. Faired curves for the SPL^ and the SPL+ al varia- 
tions at each microphone height can then be used to provide a reasonably 
accurate determination of the ”true" lossless spectrum for each test con- 
dition. 


Selective Placement of 
Reinforcements and Cancellations 

As an alternative to the computer method to correct ground reflection 
effects, use can be made of selective placement of reinforcements and can- 
cellations. This can be achieved most simply by changing the microphone 
height. For a given rig geometry having a specific source height and at a 
given radius, moving the microphone close to the ground moves the location 
of the first cancellation to higher frequencies, while moving the micro- 
phone far from the ground moves the location of the first cancellation to 
lower frequencies. Calculations to select the most desirable microphones 
heights can be made using the equations presented herein to simulate mea- 
sured SPL spectra for different test geometries. 

Another way to correct ground reflection effects is to combine near- 
er on-the-ground microphone spectra with high-above-the-ground microphone 
spectra. However, while the reflectance for measurements made by near- or 
on-the-ground microphone is constant over a large frequency range (eq. (13)), 
the value of the constant level still depends on the magnitude of the re- 
flected signal transfer function, Q. Also, while a microphone high above 
the ground produces a measured spectrum level which has its first cancella- 
tion at low frequencies, the average reflectance at high frequencies still 
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depends on the magnitude of the reflected signal transfer function* Hence 
while near- or on-ttie-ground microphone spectra can be combined with high- 
above-the-ground microphone spectra to obtain data free of cancellations 
and reinforcements, the effect of ground reflections may still be present. 
Some estimate of the value of Q should be made and the proper correction 
should be applied to both ground microphone and high-above-the-ground 
microphone spectra data before they are combined. The computer method 
discussed herein can be used to obtain the necessary value of Q. 

From the data presented herein, it is speculated that if the lossless 
SPL spectra accuracy needs to be only ±2 dB, it would be sufficient to 
assume the reflectance for a microphone near or on the ground can be ex- 
pressed by 

| T | 2 = 10 log 10 ( 1 + 1/z) 2 - 1 (17) 

and the reflectance for a microphone high above the ground can be expressed 
by 

' 1 T | 2 = 10 log 10 (1 + 1/z) 2 -1 (18) 

These assumptions could be checked using the method described herein. 

In general the frequencies of the cancellations and reinforcements 
are fixed by the facility geometry. The magnitude of the problem caused 
by reflection effects in any test will depend on the frequency range of 
interest of the scale model being tested. For studies where the frequency 
range of interest is above the location of the cancellations and reinforce- 
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merits and where only relative level changes are important, the ground 
reflection effects may not be pertinent, 

CONCLUDING REMARKS 

This study has shown that the effects of ground reflections on jet 
noise spectra can be satisfactorily corrected by a computer method. The 
computer method was based on an analysis of ground reflection effects 
using fhe concept of a reflected signal transfer function which represents 
the effects of both the reflecting surface and the atmosphere on the 
propagation of the reflected signal in a mathematical model. The com- 
puter method, besides being useful in correcting measured data, can be 
employed to predict the results of using a particular microphone array 
or test site treatment. 

Ground effects can also be avoided or minimized by appropriate 
placement of the microphone array to produce selective placement of rein-^ 
forcemeats and cancellations. The computer method can also be helpful 
in this approach. 
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APPENDIX A 
SYMBOLS 

atmospheric attenuation level, dB 

cost function, sum of square of error between measured 
and calculated sound spectra 

PSD + function level at zero Strouhal number, (eq. (B-3)), dB 

penalty function 

reflection coefficient 

total cost function 

cancellation frequencies, Hz 

velocity of sound in air, m/sec 

velocity of sound in surface, m/sec 

effective diameter, m 

Strouhal response function, (eq, (B-4)) 
frequency, Hz 

nominal one-third-octave band frequency, Hz 
reflected signal transfer function 
microphone height, m 
height of n^ source, m 

number of E terms summed to form PSD + 
number of independent noise sources 
number of E terms subtracted to form PSD + 
source index 
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OASPL 

OASPL + 

r 

PSD 

PSD + 


PSD 


+ 

cal 


Pr 

Q(s) 

r 2 »®- 3»®-4 

r 

SPL 

SPL + 


SPL 


cal 

SPI 4l 

spp m 

SP < 

St 

St 


overall sound pressure level relative to p r , dB 
lossless OASPL parametric value, dB 

normalized mean-square pressure spectral density level, dB 
PSD without atmospheric absorption and ground reflection 
effects, dB 

PSD + calculated by an equation for a real random process 
spectrum density, dB 

reference pressure 2x10”^ Pascal (2x10”^ microbar) 
magnitude of reflected signal transfer function 
reinforcement frequency 
distance from source to microphone 

one-third-octave sound pressure level referenced to p r , dB 
SPL that would be measured without atmospheric attenuation 
and in the absence of a reflecting surface (lossless 
SPL), dB 

SPL calculated from equation (16) , dB 
lossless SPL calculated from equation (B-l) , dB 
SPL measured, dB 

lossless SPL calculated from data by equation (12), dB 
Strouhal number, f(D e /V e ) 
critical Strouhal number 
j2TTf, jOO 

distance from source to microphone along direct signal 


path, m 
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s 2 

I T(f c ) | 2 

t 

u 

v e 

21 

22 


z 


Af c 
<5 (s) 

n p 


0 


2 


pc 0 

T 


distance from source to microphone along reflected signal 
path, m 

one-third-octave reflectance 
time, sec 

Strouhal number ratio, dimensionless 
effective velocity, m/sec 

air impedance, acoustical ohm-m (mks rayl) 
complex ground surface impedance, acoustical ohm-m^ 

(mks rayl) 

mean ratio of path length reflected signal traveled to 
path length direct signal traveled for all n sources, 
s 2 /si 

ground absorption coefficient 

bandwidth at one-third-octave band center frequency, 
f c (0.23 lb f c ), Hz 

reflected signal transfer function phase factor, radians 
Strouhal response excitation parameter 
angle between incident acoustic signal and normal to 
surface, radians 
see equation (18) 
density of air, kg/m^ 

acoustical impedance of air, 416 acoustical ohm-m^j (mks rayl) 
time delay, sec 

acoustic delay time from source to microphone, sec 
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0) 

Subscripts: 

cal 

i 

k 

m 

o 

P 

pen 

R 

r 

s 

T 

Superscripts 

+ 


angular velocity; 2’irf, rad/sec 
calculated 

index of numerator terms in equation (B-3) 

index of denominator terms in equation (B-3) 

measured 

in air 

parametric 

penalty 

reflection 

reference 

source 

total 

evaluated with reflection and atmospheric attenuation 
effects removed 
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APPENDIX B 

ANALYTICAL LOSSLESS SPL SPECTRA 
The lossless third octave SPL spectra equation used herein is 

OASPLp / 10 D 

S p L ca i(J.) = 10 log 10 10 Af c — + PSD cal (B-l) 


where 


OASPL+ = 10 log Y. 10 (B-2) 

P U t 10 cal 

(All symbols are defined in the symbol list.) The quantity denoted as 
PSD ca j_ is the normalized mean square pressure spectral density level. 

The following semi-empirical formula presented in Appendix E of 
reference 3 is used herein to represent the PSD^^ function. 

M a N b 

PSD cal = C o + Z_ E[u(i), n («] -£ E [udO, n p (k)] (B-3) 

l—l 1c 1 


where E is the Strouhal response function which represents a real random 
process spectral density 


E(u, r| p ) = 10 log 1Q 


n 2 + (u-i) 2 
p 

n 2 + (u+i) 

. + 1 . 

ri 2 + u 2 +l 


(B-4) 


V 
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u is the Strouhal number ratio 


u = St/Stp 


p is the Strouhal response excitation parameter, and St p is the 
critical Strouhal number. The normalization constant, C G from equation 
(B-l) and (B-2) is adjusted so that 


z 


Af a) d 
c e 

V 


PSD tal W/1 ° 


10 


1 


(B-5) 
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APPENDIX C 

THE REFLECTED SIGNAL TRANSFER FUNCTION 

The relatively large range of values of the magnitude of the reflected 
signal transfer function, Q, used in analyzing the data shown in figure 5 
can be explained by two different interpretations of the reflected signal 
transfer function (G(f c ) = Q e *^). Briefly, the reflected signal transfer 
function can be considered as providing information about the reflecting 
surface and/or as providing information about the propagation of the re- 
flected signal. The results will be interpreted in terms of these two 
factors. 

Propagation processes . - In terms of the propagation process, the 
value of the magnitude of Q indicates the amount of correlation between 
the direct and reflected signal. If the direct and reflected signals were 
not correlated, the measured spectrum would represent the sum of the squares 
of the direct and reflected signals and no cancellations or reinforcements 
would occur. 

The occurrence of a small cancellation and reinforcement effect makes 
interpretation of the SPL^ data difficult. As an example, the data taken 
at a 7.62 m (25 ft) radius for microphone heights of 1.52 and 3.89 m (5 and 
12 3/4 ft) above the asphalt surface (shown in fig. 4(a)) was duplicated by 
the computer procedure using a Q value of 0,3 (fig* 5(a)). This value of 
Q indicates the presence of cancellations and reinforcements of reduced 
magnitude. However, the data taken at 7.62 m (25 ft) radius for microphone 
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heights of 0.30 and .91 in (1 and 3 ft) above the asphalt surface (shown 
in fig. 4(a)) was duplicated using Q values of 0.5 and 0.614, Since the 
asphalt surface was the same in both cases, this indicates that the corre- 
lation between the direct and reflected signals could have been modified 
by other propagation processes than ground reflection effects. Some propa- 
gation processes which might be a factor are scattering, turbulence in the 
atmosphere and fluctuations of wing and temperature gradients with time. 
Investigation of these factors is beyond the scope of this report. 

Surface impedance . - If it is assumed that the direct and reflected 
signal correlation are unaffected by atmospheric effects, then the asphalt 
surface impedance can be a factor. The computer in correcting the data 
calculates the reflected signal transfer function, G = Q e . Figure 8 
shows a typical plot of the complex number G for a given microphone 
distance. For all of the distances covered, the magnitude of G ranged 
from 0.3 to 0.8 and the average value was near 0.6. 

From equation (7) the ratio of the impedance of the asphalt to that 
of air is 



1 + Q 16 \ 
e 

cos 02 

1 + Q e 16 1 

COS 0J_ 


(Cl) 


Using the values shown in figure 8 for G and letting cos ©2 - 1 the 
value of was found. Typical values of ^/^l are P^ otte ^ 

figure 9. For all the cases studied the magnitude of ran 8 ed fro™ 

10 to 25 for microphones at 15.24 and 22.66 m (50 and 75 ft) and the aver- 
age value was about 17 . 



31 


Reflectance . - Figure 10 shows the reflectance calculated using the 
assumption Q = 1 (perfect reflection) and using the computer selected 
value of Q. Comparison indicates that the cancellations and reinforce- 
ments are larger when Q = 1. The value of Q selected by the computer 
decreases the magnitude of the cancellations and reinforcements and also 
introduces a phase shift. From equation (3), the first cancellation 
occurs when 

cos (27rfT n - 6) = -1 

or for 6 = 0, at 

f - JL. c ° 

2T n 2(s2~S]_) 

When 6 is near zero, the first cancellation occurs at 

f _ 6 +tt C o 
r c ~ 27 t Ts -s^)" 

Thus the computer as shown in figure 10 selected small positive values of 
6 for microphone heights of 0.91 and 3.89 m (3 and 12 3/4 ft). This re- 
sulted in the shift of the first reinforcement about one third octave band 
higher than when 6 was zero. This phase shift may be due to the impedance 
of the asphalt or to changes in path length from temperature or wind gradi- 
ents. Thus the impedances shown in figure 9 are probably in error. 
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TABLE X. - MICROPHONE LOCATIONS 


Microphone 

Microphone 

Microphone 

Angle, 

Radius , 

Height, 

Degrees, 

r 

Jkx 


m ft 

m ft 



7.62 m (25 

ft) from 

nozzle 

exit 


160 

6.55 

21.51 


0.0061 

0.02 

155 

6.71 

22 


0,3048 

1 

150 

7.01 

23 


0.9144 

3 

145 

7.32 

24 


1.5240 

5 

145 

7.62 

25 


3.89 

12.75 


15.24 m (50 

ft) from 

nozzle 

exit 


150 

14. 63 

48 


0.0061 

0.02 

150 

14.63 

48 


0.3048 

1 

145 

15,24 

50 


0.9144 

3 

145 

15.24 

50 


1.5240 

5 

145 

15.24 

50 


3.89 

12.75 


22.86 m (75 

ft) from 

nozzle 

exit 


145 

22.86 

75 


0.0061 

0.02 

145 

22.86 

75 


0.9144 

3 

145 

22.86 

75 


1.5240 

5 

145 

22.86 

75 


3.89 

12.75 


7.62 m (25 

ft) from 

nozzle 

exit 


125 

7.62 

25 


3.89 

12.75 


100 7.62 25 3.89 12.75 
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Figure 1. - Schematic of nozzle and microphone geometry relative 
to asphalt surface. 
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Figure 2. - Diagram of computer method. 

















1C) LOSSLESS SPECTRA 


Figure 3. - Example of measured and computer-calculated 
spectral results. Microphone height, 0. 3048 m; source 
height, 3. 89 m; radius, 7. 62 m. Nozzle exhaust veloc- 
ity 301 m/sec; directivity angle, 145°. Q = 0. 605, 

6 = 6. 36. 
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Figure! - Concluded, 
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Figure 4, - Effect of microphone height on measured spec- 
trum for various microphone radial distances. Directiv- 
ity angle, 145°. 


MICROPHONE HEIGHT, 



Figure 4 - Concluded. 
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(a> MICROPHONE DISTANCE, 7.62 m. 




MICROPHONE HEIGHT, 

Q 

M 2 


h 0 . m 


dB 

0 

0.0061 

1.0 

6 

□ 

0.0061 

0.637 

4.3 


0 91 

0.614 



1. 52 

0.664 


— 

3.89 

0.534 




10 MICROPHONE DISTANCE, 22.86 m. 


Figure 5. -Comparison of lossless sound pressure level 
spectra based on corrected measured data taken at 
ground level and at various microphone heights for 
several microphone radial distances. 
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lb) h 0 - 1. 52 mj Q - 0. 664 

Figure 6. - Comparison of measured and lossless sound 
pressure level spectra based on correcting measured 
data for various microphone heights. Microphone 
radial distance, 22 . 56. Exhaust velocity, 301 m/sec. 



Figure 6. - Concluded. 


66Z8-3 




Figure 7. - Comparison of lossless sound pressure level 
based on measured data, SPL^,, and os calculated 
from curve fit, SPl™ L . Microphone radial distance. 
22.86 m. Exhaust velocity, 301 m /sec. 
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Figure 8. - Reflected signal transfer 
function, G(t c ) for various micro* 
phone heights. Microphone radi- 
al distance, 22. 66 m. 
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Figure 9. - Surface impedance relative 
to air, for various heights. 

Microphone radial distance, 22. 86 m. 
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Figure 10. - Comparison of reflectance calculated for per- 
fectly reflecting ground (Q ■ 1) with reflectance calcu- 
lated by computer selected constants. Microphone 
radial distance from nozzle, 22, 86 m. 
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